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We study the structure, the formation and binding energies and the transfer levels of the zinc- 
phosphorus vacancy complex [Znj n -Vp] in Zn doped p-type InP, as a function of the charge, using 
plane wave ab initio DFT-LDA calculations in a 64 atom supercell. We find a binding energy of 
0.39 eV for the complex, which is neutral in p-type material, the 0/-1 transfer level lying 0.50 
eV above the valence band edge, all in agreement with recent positron annihilation experiments. 
This indicates that, whilst the formation of phosphorus vacancies (Vp 1 ) may be involved in carrier 
compensation in heavily Zn doped material, the formation of Zn-vacancy complexes is not. 

Regarding the structure: for charge states Q — +6 — » —4 the Zn atom is in an sp 2 bonded 
DX position and electrons added/removed go to/come from the remaining dangling bonds on the 
triangle of In atoms. This reduces the effective vacancy volume monatonically as electrons are added 
to the complex, also in agreement with experiment. The reduction occurs through a combination of 
increased In-In bonding and increased Zn-In electrostatic attraction. In addition, for certain charge 
states we find complex Jahn- Teller behaviour in which up to three different structures, (with the In 
triangle dimerised, antidimerised or symmetric) are stable and are close to degenerate. We are able 
to predict and successfully explain the structural behaviour of this complex using a simple tight 
binding model. 

PACS numbers: 71.55.Eq 61.72.Bb 71.70.Ej 71.15.Dx 



I. INTRODUCTION 

Zn doped p-type InP is one of the most common ma- 
terials in use within optoelectronics. The Zn sits substi- 
tutionally within the In sublattice (Znj n ) where it has a 
shallow acceptor level. A well known limitation in the 
use of Zn as a p-dopant in InP is the saturation of the 
hole concentration in the mid 10 18 cm~ 3 range. Above 
this further increases in Zn concentration do not trans- 
late into increases in hole concentration. There are sev- 
eral suspected causesi for this, in particular: a) increases 
in the concentration of interstitial zinc (Zn;), b) phase 
separation and c) the formation of other compensating 
defects, especially phosphorus vacancies 1 (Vp) and com- 
plexes of zinc with Vp . In this paper we will examine the 
properties of Zn-Vp complexes, using ab initio methods. 

The formation of one such Zn-Vp complex has recently 
been studied^ by Slotte et al. in Zn doped Czochralski 
grown samples which have been post-growth annealed in 
the temperature range 300-600 K. Using positron annihi- 
lation they find that the complex has at least two forms. 
The most stable at room temperature^ has an effective 
vacancy volume larger— than that of a free Vp . A bind- 
ing energy on the order of 0.1 — > 0.4 eV is anticipated™ 
(prior to detailed measurements). When the measure- 
ment temperature is raised they find that the complex 
undergoes a transition to a form with a smaller vacancy 
volume, similar to that of the free vacancy. It has been 
suggesed that the more stable form is neutral^ and has 
a DX structure^ which becomes non-DX when excited to 
the -1 charge state. They findA the 0/-1 transfer level 
0.2±0.1 eV above the valence band edge. 



A number of theoretical studies of the free phos- 
phorus vacancy have been performed^ii using Density 
Functional Theory (DFT). The free vacancy has been 
shown§i£ to be a strong Jahn- Teller defect, indeed, to be 
a "negative U" centre in which the neutral charge state 
Vp is thermodynamically unstable relative to the +1 and 
-1 states. Hence a direct change Vp 1 — » Vp 1 occurs as a 
function of increasing electron chemical potential (Fermi 
level). This is due to strong lattice relaxation in the -1 
charge state. The vacancy is surrounded by a tetrahedron 
of four In atoms, each with a (partially filled) inwardly 
pointing dangling bond. In the -1 charge state a pair of 
dimers is formed: two opposing sides of the tetrahedron 
are reduced in length relative to the remaining four. The 
energy gained by this overcomes the Coulomb repulsion 
between the additional electrons, allowing two electrons 
to be transfered onto the vacancy simultaneously. 

In this paper we present what is, to our knowledge, 
the first ab initio^ ab initio theoretical study of the 
[Znj n -Vp] complex, which we propose to be that ob- 
served in the positron annihilation experiments. We will 
first describe the method in section [D] We will then 
report detailed results for the geometry and vacancy vol- 
ume of the complex, both for the ground state structures 
and for various metastable structures in section llTll fol- 
lowed by a description of how these vary with changes in 
the charge state of the complex in section ITVl In section 
0we will also propose a simple tight binding model which 
is able to predict almost all of the structural properties 
which wc find. Finally, in section IVII we will conclude, 
relating our results to those from position annihilation 
and the issue of Zn compensation. 
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TABLE I: Comparison of the k-point convergence of the elec- 
tronic and structural contributions to the formation energy, 
Ed- Convergence of the electronic contributions is given 
by the change, A c/ii (N-M), in the unrelaxed formation en- 
ergy when increasing the k-point grid from NxNxN to and 
MxMxM. Convergence of the structural contributions is given 
by A fl (N-M), the equivalent change in relaxation energy 
en(N). See main text for full definitions. 
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II. 


CALCULATIONAL DETAILS. 





We use planewave ab initio DFT— within the Local 
Density Approximation (LDA) together with ultrasoft 
pseudopotential o 9 ! 10 . Calculations are performed using 
the VASP codeii. Before doing calculations involving 
defects we first optimise the lattice parameter subject to 
LDA. We find a value of 5.827 A, compared to 5.869 A ex- 
perimentally, giving us a band gap of 0.667 eV, compared 
to 1.344 eV in experiment. 

For defect calculations we use a 64 atom simple cubic 
supercell and allow all atoms not located on the surface 
of the cell to relax. This restriction is included to trun- 
cate the (spurious) elastic interactions between adjacent 
supercells. No restrictions are placed upon the symme- 
try of the relaxations. When calculating total energies 
of supercells containing charged defects we also include a 
uniform neutralising background charge. The key quan- 
tity is the formation energy of the defect, Ed, defined^ 



as 

E d = £ T (defect Q ) - £ T (bulk) + ^ mm + Q (e v + e/) 

i 

(i) 

where E T (defect^) and i? T (bulk) are the total energy 
of the InP supercell with and without the charge Q 
defect. Both are calculated with the same values of 
planewave cutoff, k-point grid, etc, in order to make 
use of cancellation of errors. The defect is formed by 
adding(removing) m atoms, each with chemical poten- 
tial fj,i, and by adding (— Q) electrons, whose chemical 
potential e/ is measured from the valence band edge e v . 
We then adjust the various calculational parameters to 
converge E&. We find that it is sufficient to use pseudopo- 
tentials in which Zn d-electrons are treated as valence 
electrons but In d-electrons are left in the core. We find 
that a planewave cutoff energy of 200 eV is enough to 
converge Ed to 0(0.01 eV). Real space projection opera- 
tors are used, with the cutoff again chosen to give errors 
below 0(0.01 eV). 

The k-point convergence is more interesting, since we 
find different convergence behaviour for the electronic 
and structural contributions to the formation energy 
Ed(K), where we use a KxKxK Monkhorst-Paclsi^ k- 
point grid. The error in the electronic contribution to 
Ed(K) is very close to the error in the unrelaxed forma- 
tion energy, E^ R (K), which we can estimate as 

A UR (N-M) = E% R (N) - E% R (M) (2) 

(We shall use {N,M} = {2,4} and {4,6} to estimate the 
errors in E^ R (2) and E^ R (A) respectively.) We cannot, 
of course, calculate the structural contribution directly, 
so we first define the relaxation energy 

e R (K) = E R (K)-EV R (K) (3) 

where E R (K) is the relaxed formation energy. We then 
estimate the convergence of c_r(K) using 

A*(N-M) = e B (N) - e H (M) (4) 

In table |T| we list values of A*™ (2-4), A ra (4-6) and 
A R (2-4) for a range of charge states of various differ- 
ent defects in InP. We find that when atomic relaxations 
are ignored, a 2x2x2 k-point grid produces errors below 
0(0.01 eV), provided Q = 0. (Column A UR (2-4) of ta- 
ble □) When Q/0, however, we find that the error 
introduced by truncating at 2x2x2 increases with \Q\. 
Instead, a 4x4x4 grid is required: column A UR (4^6) in- 
dicates that this is sufficiently converged. On the other 
hand, the relaxation energy £ij(K) converges faster, see 
column A R (2-4) in the table. Indeed, £r(K) is converged 
for a 2x2x2 Monkhorst-Pack grid, for any value of Q. In 
other words the structural contribution to E d appears to 
converge faster with k-point grid than the electronic con- 
tributions. This is largely due to a double cancellation 
of errors when we calculate sr from values of Ed which 
already contain a cancellation of errors. It means that, 
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TABLE II: Transfer levels of [Znj n -Vp] , in the 64 atom cell, 
relative to the valence band edge. 



Level 


Energy 


+6/+5 


-0.51 eV 


+5/+4 


-0.51 eV 


+4/+3 


-0.47 eV 


+3/+2 


-0.44 eV 




-U.ol ev 


+1/0 


-0.32 eV 


0/-1 


0.50 eV 


-1/-2 


0.70 eV 


-2/-3 


0.95 eV 


-3/-4 


1.05 eV 


-4/-5 


1.19 eV 


-5/-6 


1.18 eV 



TABLE III: Geometry of the most stable form of [Znj n -Vp] 
as a function of Q. Distances in A, volumes in A 3 . Most struc- 
tures (Str.) found are symmetric (SY) (all In-In distances 
equal, all Zn-In distances equal.) Where the In-In (and like- 
wise Zn-In) distances differ it is found that two remain equal 
(In-In(2)) and one is different (In-In(l)) by amount A(In). 
(In — In(l)) < (In — Jn(2)) corresponds to the formation of a 
dimer (DM) and (In — Jn(l)) > (In — In(2)) to an antidimer 



(AD). 


Q 
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Volume 


Str. 


A(In) 






(1) (2) 


(1) (2) 
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2.29 


4.97 


4.95 


14.36 


SY 




+5 


2.29 


4.97 


4.96 


14.42 


SY 




+4 


2.29 


4.98 


4.97 


14.49 


SY 




+3 


2.29 


4.97 


4.99 


14.53 


SY 




+2 


2.29 


4.99 


4.97 


14.57 


SY 




+1 


2.30 


4.69 


4.36 


10.86 


SY 




+0 


2.30 


4.41 


3.69 


7.58 


SY 




-1 


2.30 


4.23 4.31 


3.64 3.37 


6.52 


AD 


7.9% 


-2 


2.31 


4.13 4.20 


3.52 3.19 


5.81 


AD 


10.3% 


-3 


2.32 


3.95 4.10 


2.99 3.28 


5.23 


DM 


8.9% 


-4 


2.32 


3.96 


3.08 


4.85 


SY 




-5 


2.43 


3.21 


3.04 


3.59 


SY 




-6 


2.51 


2.99 


2.95 


3.09 


SY 





although we should use a 4x4x4 grid to obtain an accu- 
rate value for E R , we can save much calculation time by 
estimating E R (A) as 

££(4) « EV R (4) + e R (2) 

= E u d R (A) + E R (2)-E u d R (2) (5) 

The resulting values of E R (A) are converged to at least 
0.02-0.04 eV, and usually to an order of magnitude or so 
better. Hence we do most relaxations with a 2x2x2 k- 
point grid, only checking key ones at 4x4x4. Charge and 
Electron Localisation Function M (ELF) plots are, on the 
other hand, calculated using a 4x4x4 k-point grid at the 
2x2x2 related structures. 



TABLE IV: Geometry of structural excitations of 
[Znj n -Vpl ^ as a function of Q. Energies are given 
relative to the energy of the most stable structures, as listed 
in table ITTT1 See table UTTI also for meaning of other collumns 



and notation. 


Q Energy Zn-P 


Zn-In 


In-In 


Volume 


Str. 


A(Jn) 


[eV] 




(1) (2) 


(1) (2) 








+4 1.00 


2.30 


4.44 


3.74 


7.85 


SY 




+3 0.97 


2.30 


4.44 


3.76 


7.91 


SY 




+2 0.93 


2.30 


4.44 


3.72 


7.76 


SY 




+1 0.12 


2.30 


4.42 


3.71 


7.67 


SY 




-1 0.011 


2.30 


4.29 4.27 


3.41 3.50 


6.57 


DM 


2.7% 


-1 0.012 


2.30 


4.28 


3.47 


6.58 


SY 




-2 0.03 


2.31 


4.18 4.17 


3.22 3.35 


5.85 


DM 


3.7% 


-3 0.02 


2.31 


4.05 4.10 


3.41 3.06 


5.26 


AD 


11.6% 


-3 0.05 


2.31 


4.08 


3.18 


5.31 


SY 




-5 0.02 


2.34 


3.67 


3.05 


4.32 


SY 




-6 0.19 


2.34 


3.63 


3.05 


4.26 


SY 





III. THE NEUTRAL [Zn In -V p ] COMPLEX. 
A. Formation and Binding Energies. 

In order to calculate absolute values for the formation 
and binding energies we need to know the chemical po- 
tentials A*j n , Hp and Mj n p- Fully converged calculations 
for In, P (black phosphorus) and InP (the latter in its 
basic 2 atom FCC cell) give the bulk values 20 -3.2 70 eV, 
-6.028 eV and -9.728 eV respectively. For pure InP we 
have, at thermodynamic equilibrium, H\nP = HTn + Hp- 
In In rich material we then have /xj n = -3.270 eV leading 
to Hp = -6.459 eV, whilst in P rich material ^p = -6.028 
eV and hence h\u = -3.700 eV. We now use the values 
/xj n = -3.485 eV and Hp = -6.243 eV, corresponding to 
exactly stoichiometric conditions. Properly, doping with 
Zn should alter these h values, but even a Zn concentra- 
tion of 5.10 18 cm~ 3 only corresponds to replacing 1 In 
in 2000 with Zn. We thus ignore the effect, and use the 
fully converged bulk value -1.891 eV for HZn- 

Since we are primarily interested in strongly p-type 
material, we assume that the electron chemical potential 
lies at the valence band edge, meaning that e/ = 0. Using 
fully relaxed calculations in the 64 atom simple cubic cell, 
we find that the most stable charge states of the Znj n and 
Vp are -1 and +1 respectively, whilst for [Znj n -Vp] it 
is 0. The formation energies are Ed(7m^ ) = 0.32 eV, 

^(V^ 1 ) = 1.96 eV and E d ( [Zn In -Vp] °) = 1.89 eV. For 

Vp we find a free vacancy volume of 6.5 A 3 , where we 
calculate only the volume of the tetrahedron formed by 
the four In surrounding the vacancy. Vp 1 is invisible to 

positrons, so we will also report here the volume of Vp 1 , 
which is 4.6 A 3 . More detailed results about Vp and 
Znj n will appear elsewhere^ 

We can evaluate the binding energy, E^, of the complex 
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as 

£ b ([Zn In -Vp]°) = E^Zn^+E.iV+^-E^Zn^-Vp] ) 

(6) 

This gives us a value of 0.39 eV, which agrees with the 
0.1 — * 0.4 eV estimate from the positron annihilation 
experiments^. 

B. Structure. 

In figure we show the relaxed structure of the 
[Znj n -Vp] complex in the Q = charge state. The Znj n 
is bonded to three P atoms only and has relaxed back into 
a DX-like position in the same plane as those P atoms. 
The resulting Zn-P distances are 2.31 A, which is almost 
identical to the 2.34 A value we find for the shorter Zn-P 
distances in a similar DFT-LDA calculation for Z113P2 
and various other Zn-P binary compounds. -~ The three 
In atoms on the other side of the Vp have moved together 
slightly, giving three equal In-In distances of 3.69 A, com- 
pared to an ideal LDA bulk next nearest neighbour dis- 
tance of 4.12 A. The Zn-In distance is rather longer than 
the bulk value, however, at 4.41 A. This gives a volume 
for the vacancy of 7.58 A 3 , where we calculate simply the 
volume of the tetradehron formed by the three In and the 
Zn. In the remaining figures most of the atoms in the cell 
are omitted for clarity, and lines are drawn around the 
edges of the tetrahetrdon to guide the eye. In the case of 
Zn-In lines they do not represent bonds. This can be seen 
clearly by examining the ELF— plots in fi gure[3 There 
is no ELF density at all between the Zn and the In tri- 
angle, indicating that there is no bonding between them. 
Instead, the Zn has three sp 2 hybrid orbitals bonded to 
the surrounding P, these bonds being filled. There is no 
density related to the Zn at all along the line perpendic- 
ular to the Zn-P plane (the (111) axis). This indicates 
that the single Zn 2p z orbital left over from the formation 
of the sp 2 hybrids is empty. 

IV. THE STRUCTURE OF [Zn In -V p ] AS A 
FUNCTION OF CHARGE STATE. 

A. Bond lengths and vacancy volumes. 

In table^we give the positions of a number of transfer 
levels of the complex. The structures of all of the states 
we have examined (+6 — > -6) are very similar, as listed 
in table EH For all Q > -5 the Zn lies in the DX- 
like position. The three In atoms remain unbonded to 
Zn, but both the In-In distances and the Zn-In distances 
shorten with increasing negative charge. This reduces 
the vacancy volume for each additional electron. This 
holds until Q = — 5, when the Zn 2p z orbital starts to 
be filled, leading at Q = —6 to a non-DX position for 
the Zn atom, with the Zn-In distance only 1.6% longer 
than the In-In. The transfer level into this state lies 



1.18 eV above the valence band edge. With the LDA 
band gap being only 0.67 eV this places the level about 
0.51 eV above the conduction band edge. This orbital 
also displays a negative U effect: the Q = — 5 structure 

- which lies half way between DX and non-DX forms - 
is thermodynamically unstable, so the transfer level -5/- 
6 lies just below the -4/-5 level. Hence we anticipate 
instead a -4/-6 level about 0.5 eV above the conduction 
band edge. 

In addition to this, for Q = — 1 — > —3, Jahn- Teller 
effects are apparent, in which one In-In distance is short- 
ened (dimerised) or lengthened (antidimerised) relative 
to the other two. These remaining In-In distances always 
remain equal to at least O.OOlA. 

In table HVl we show the various metastable structures 
(structural excitations) which we find. (An extensive 
search suggests that there are no other metastable struc- 
tures besides these— The fact that the Zn 2p z level lies 
well inside the conduction band is further underlined by 
the metastable structures seen for Q = — 5 and Q — — 6. 
For both of these we find a DX-like metastable structure 
very similar to the ground state structure of the Q = —4 
state. This indicates that one or two electrons (respec- 
tively) have been transfered from the complex into the 
conduction band, leaving the complex with an empty 
Zn 2p z and additional charge elsewhere in the super- 
cell. The energy splitting between these DX and non-DX 
forms is therefore expected to change significantly (in- 
deed, to switch over) in larger supercells where less de- 
fect level dispersion is seen. A similar effect is seen in the 
Q = +1 — > +4 charge states, where a metastable struc- 
ture is found which closely resembles the Q — structure. 
Since the hole levels corresponding to these charge states 
lie below the valence band edge, these structural excita- 
tions amount to exciting holes off the complex, leaving 
it in the Q = charge state, with additional charge else- 
where. (An equivalent structural excitation for Q = +5 
or +6 has not been found - presumably the excitation en- 
ergy involved is too great for even metastability within 
this supercell.) 

All of the other excited structures found in the cal- 
culations relate to variations in Jahn- Teller distortions 
and affect only the In-In and Zn-In bond lengths. They 
involve no significant changes in vacancy volume, Zn po- 
sition or electron occupancy of the complex (see next 
subsection). The energies involved are very small indeed 

- on or below the order of room temperature thermal 
energies. Room temperature experiments are therefore 
likely to see averages over these structures. Three pos- 
sible structures are found to be stable - one In-In bond 
dimerised (the other two equal), one antidimerised or all 
equal. How the symmetric structure comes to be stable 
despite the Jahn- Teller effect will be returned to in sec- 
tion All three structures are found to be stable for 
Q = — 1 and Q = —3, but for Q = —2 the symmetric 
structure is not seen. For Q = and +1 we find only the 
symmetric structure. 

The energy separations between these structures are 
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TABLE V: Ab initio results for the Lowest Lying structures 
(marked L.L.) of [Zn In -Vp] Q for Q = +2 -» -4, together 
with structural excitation energies in electron volts. DM = 



dimerised, SY = symmetric, AD = antidimerised. "-" indi- 
cates that the structure is not stable. 

Q DM SY AD 

+2 - L.L. 

+1 - L.L. 

- L.L. 

-1 L.L. 0.0004 0.0001 

-2 0.0002 - L.L. 

-3 0.0097 0.0104 L.L. 

-4 - L.L. 



on the limits of what we expect to resolve using only the 
2x2x2 Monkhorst-Pack k-point grid together with equa- 
tion El As a result we have re-calculated the energies 
and structures of all of these states using a 4x4x4 k-point 
Monkhorst-Pack grid, and, for some, even a 6x6x6 grid. 
We find that all of the structures reported in table IIIII 
andflvl remain stable, although the degree of dimerisa- 
tion or antidimerisation changes in some cases. The order 
of stability for the three structures changes for Q = — 1 
and Q = —3, however, and in all charge states they move 
closer to degeneracy. These changes are illustrated in 
table E| 



B. Charge density differences. 

In figureQIwe show the instantaneous change in charge 
density when an electron is added to (upper row) or re- 
moved from (lower row) the Q = charge state. In other 
words, we first take the Q = relaxed positions and cal- 
culate with the Q — ~ 1 or Q = +1 total charge. We 
then subtract from this new charge density the original 
charge density from the relaxed Q — calculation. The 
isocharge surfaces plotted are for 80%, 50% and 20% of 
the peak value. It is clear that an electron added goes to 
a localised orbital on or near the In triangle. An electron 
removed, however, comes from a delocalised state. We 
note, however, that even for the 20% peak value isocharge 
surface for 0/-1, the total charge enclosed corresponds to 
only 0.48 electrons. Hence, even for the Q = — 1 charge 
state, most of the added electron has an apparently de- 
localised character. This is, however, a result of the use 
of a relatively small supercell: the In-In distance is 3.7 
A within the cell, but about 9.2 A between cells. This 
produces a significant fake dispersion for the defect lev- 
els, and means that the transfer levels given in table ITT1 
should not be trusted too heavily. 

In figure 0] we show the 50% peak isocharge surface 
for the instantaneous charge density change upon adding 
or removing electrons from the symmetric structures for 
Q = +1, and -1. The orbital filled when adding an 
electron to Q = and that emptied when removing an 
electron from Q = — 1 are essentially the same. Clearly 



Q = — 1 is a localised state of the [Znj n -Vp] complex. 
Adding an electron to the Q = — 1 state gives a very 
similar shaped charge density - the opposite spin state of 
the same orbital. 

Removing an electron from Q = empties a delocalised 
valence band state, as can be anticipated from the trans- 
fer level being found inside the valence band in the pre- 
vious section. However, adding an electron to the (fully 
relaxed) Q = +1 fills an In triangle localised orbital, none 
the less. Clearly Q — + 1 also corresponds to a localised 
state of the complex. 

Going beyond +1 and -1, figure |U also indicates a pat- 
tern which we in fact find for all charge states from +6 
to -6: the electrons added or removed at the transfer 
levels +2/+1 and +1/0 come from essentially the same 
orbital, which we shall call <&o. The same is also try for 
the orbitals corresponding to the transfer level pairs 0/- 
1 and -1/-2, ($i), -2/-3 and -3/-4 ($ 2 ) and -4/-5 with 
-5/-6 ($3) Rather than plot all of them we show only 
one example of each orbital (transfer level pair). For the 
first three of these orbitals the shape is best seen in a top 
view, looking down the (111) axis. This is shown in figure 
[5] for the symmetric, dimerised and antidimerised struc- 
tures. A clear relationship emerges: the dimer orbital 
for 0/-1 is very similar to the antidimer orbital for -2/-3. 
Likewise, the dimer orbital for -2/-3 and the antidimer 
orbital for 0/-1 are the same. Furthermore, a symmetric 
combination (addition) of these two would give an or- 
bital very much like that of the symmetric -2/-3 orbital. 
An antisymmetric combination (subtraction) would give 
something similar to the symmetric orbital for 0/-1, al- 
though in this case the similarity is not complete. This 
may be due to differences in the degrees of dimcrisation 
or perhaps to limitations in the calculation. Despite this, 
a clear relationship between the orbitals filled/emptied at 
these transfer levels in the three nearly degenerate struc- 
tures is present and will be considered in more detail in 
section Ivl 

For the orbital associated with -4/-5 and -5/-6 we show 
a side view in figure [S] It is the Zn 2p z orbital, as ex- 
pected based upon the structures reported in the previous 
section. Confirming this now allows us to fully explain 
the volume dependence of the complex, as shown in table 
lllll The volume is reduced monatonically from Q = +2 
to Q = —6. For Q = —4 — * —6 this is due simply to the 
change from DX to non-DX - as the Zn 2p z orbital fills 
up, the sp 2 bonding to the P atoms is reorganised, with 
a normal sp 3 dangling bond developing. This can then 
bond on a more or less equal footing with the In dangling 
bonds. Then from Q = +2 — > —4 the In-In distances are 
reduced as electrons enter the In dangling bonds, increas- 
ing the bonding between them. The volume reduction, 
however, is due only partly to this: partly it is due to the 
concurrent reduction in Zn-In distances, even though no 
bonds form between the Zn and In until Q = — 5. This 
reduction is electrostatic in origin. The empty Zn 2p z 
orbital, together with the relatively high ionicity of the 
Zn-P bonds, (see figure |2|) leaves the positive Zn core 
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slightly less well screened along the direction pointing 
towards the three In atoms. The negative charge on the 
In atoms, meanwhile, increases as electrons are added, so 
the Zn-In distance is reduced electrostatically. Thus the 
reduction in effective vacancy volume at the 0/-1 transfer 
level, believed to have been observed in the positron an- 
nihilation studies^ is due to a combination of the In-In 
bonding and Zn-In electrostatic attraction, rather than 
to a DX to non-DX transition. 

Regarding the other structures in tables IIIII and IIVI 
we find that, as anticipated from the structures them- 
selves, the orbitals filled by adding further electrons 
to [Znj n -Vp] or removing further electrons from 

[Znj n -Vp] +2 are all delocalised. We show only one ex- 
ample in figure El namely removing an electron from 
+5. The orbital filled by removing an electron to the 
metastable shrunken structure of Q = +2 (also in figure 
|SJ) is identical to that associated with removing an elec- 
tron from the ground structure of Q = +0. This confirms 
the explaination given above of these shrunken structures 
in terms of holes excited off into delocalised valence band 
states. 



V. STRUCTURAL MODEL FOR THE [Zn In -V p ] 
COMPLEX. 

A. Symmetry considerations. 

Since, for Q > —5, the Zn is in a DX-like position and 
the Zn 2p z is much higher in energy than the other or- 
bitals involved, we omit it, and base our model on the 
assumption that all additional electrons added/removed 
from the complex go to/come from the triangle of sp 3 hy- 
brid dangling bonds on the In atoms. These start out life 
transforming as the + (s) and 1~ (p) irreducible repre- 
sentations of the O3 group for the free atoms, which de- 
scend to the (1 dimensional) A\ and (3 dimensional) T2 
representations of the tetragonal group Td when placed 
inside a zinc blend structured crystal. (See the level split- 
ting diagram in figure 0) The point group of a free va- 
cancy (without the Zn) is the same - Td - and so the 
bound states inside the vacancy are known to also be 
described^ by the same A\ and T 2 representations of Td 
- giving a deep symmetric level and a three-fold degener- 
ate excited level. (Jahn- Teller distortions lifting this de- 
generacy lead to the negative U effects for the free Vp.) 
In [Znj n -Vp] the Zn further lowers the point group to 
G% v . The A\ representation of Td is projected onto the 
A\ representation of Cs v , whilst the Ti representation 
splits, giving a second A\ representation and a 2 dimen- 
sional E representation. The two A\ representations can 
mix, so that the hybridisation on the Zn is reorganised, 
to give the high (unoccupied) Zn 2p z orbital observed in 
our calculations, plus a new totally symmetric level con- 
structed primarily from the three remaining In dangling 
bonds. The E representation corresponds to a two-fold 



degenerate excited state, also constructed primarily from 
In dangling bonds. In the neutral state the lowest A\ 
level is full and the others are empty. For charge states 
Q = — 1 — * — 3 this degeneracy must also be lifted, usu- 
ally by a further structural Jahn- Teller effect. This can 
happen in two ways - via the formation of one stronger 
or one weaker bond, leading to the C\ v point group, or 
via a process reducing the symmetry to C3, splitting the 
E level into / = ±1 components. The latter requires a 
rigid rotation of the In triangle around the (111) axis, al- 
ternately shortening and lengthening the In-P bonds in 
the plane perpendicular to (111). This is not likely to 
be favoured, however, since the energy gained would be 
very small. Indeed, we do not see any tendency towards 
this in our ab initio structures. 

For certain Q it is also possible for the degeneracies to 
be lifted by spin interactions. In the Q — — 1 state, where 
the In triangle contains 3 electrons, it would be possible 
to localise one electron on each In - making the lower 
Ai and E represenations degenerate - and then form an 
5 = | high spin state. This would required the In sp 3 
dangling bonds to be spatially rather small (localised) 
and hence local Coulomb interactions to be rather strong. 
This does not appear to be the case, however. In addi- 
tional calculations using the Local Spin Density Approx- 
imation (LSDA) to search for it we find no evidence of 
such a state being stable. Similarly, in Q = —2, with 
4 electrons present, an S = 1 high spin state could be 
formed in the half filled Ti. Again, we find no indica- 
tions that this is stable. 

Instead, it seems that the symmetric structures for 
Q = — 1 — > — 3 are stablised by weak Coulomb interations 
with no localisation. We can estimate the differences be- 
tween the classical spatial overlaps of the three orbitals. 
Specifically, we take the charge contained in the overlap 
regions of the isocharge surfaces for charge density dif- 
ferences, (as plotted in figures EEJ) ) renormalised to a 
whole electron to allow comparisons, and evaluate what 
fraction of the charge overlaps classically. Using the 50% 
isocharge surfaces in figure [SJ we find that the highest 
orbital (corresponding to transfer levels -2/-3 and -3/-4) 
has 1.58 times as much overlap with the lowest (trans- 
fer levels +2/+1 and +1/0) as the middle orbital (0/-1 
and -1/-2) has with the lowest. (The overlap between the 
middle and the highest is about 29 times larger again.) 
Using the 20% isocharge surface, where more of the ac- 
tual charge is accounted for, but where the spurious de- 
fect level dispersion is starting to wash out the effect, we 
find that the lowest orbital has 1.62 times as much over- 
lap with the highest as it has with the middle one. These 
numbers indicates that in the symmetric geometries the 
E level is split into orbitals having less and more overlap 
with the deep A\ level. These levels are thus split by 
Coulomb repulsion from the electrons filling the deep A\ 
level. 



7 



B. A tight binding model. 

The Jahn- Teller behaviour found in the ab initio stud- 
ies can be explained using a simple tight binding model. 
To do this we need only to consider the leading order 
effects of changes in the In-In distances. We make no 
assumptions yet about the form of the degenerate associ- 
ated orbitals <f>i and we consider only the "hopping" over- 
lap matrix elements tij which move electrons between the 
orbitals on In^ and on In^. The Hamiltonian is simply 

where cj a (ci. a ) creates (annihilates) a spin a electron at 
i, the sums being over the three In atoms. We make the 
selection (23 = t + S and tyi = t — | — ti 3 . Hence 5 > 
corresponds to forming a single In2-In3 dimer (plus two 
antidimers) and S < to a single In2-In3 antidimer (plus 
two dimers). (The — | terms balancing the +6 allow di- 
rect comparisons between 5 > 0, = and < 0.) We omit 
Coulomb repulsion (we assume that U <C t) and also 
the spin-spin interactions which would lead to high spin 
states. On the other hand, we will treat the symmet- 
ric structures as (potentially) metastable, thus implicitly 
assuming the presense of (at least a) vanishing Coulomb 
repulsion between the E levels and the deep Ai level. 
More importantly we also omit elastic effects from the 
lattice. To leading order, 5 is linear in the change in 
In-In distance, but the elastic energy is quadratic. We 
can therefore write the elastic energy for an individual 
bond between In* and In^ as XSf^t, where A is (related 
to) an elastic constant. Hence for both the dimer and 
antidimer states described above the elastic energy is 
« XS 2 t + 2A (| ) 2 t = 2Mfl. However, we omit this term 
at present, solving only equation 

In the dimer case (6 > 0) the solution is: 



Bk= + ^ + £ (8) 

E 2 = 1 + 6 
in units of t, (expansions to 0(S 2 ) only) where 

F{6) = V3 V3 -28 + S 2 
5 2 

«3-<5+ y (9) 

We write the wavefunctions in the form $ = (01,02,0(3) 
where {a^} are the coefficients of {4>i} in <&. The wave- 
functions associated with Eq, E\ and E2 are thus 

*x= (i±±±^,l,l) „(-2-J-^,l,l) (10) 
$2= (0,1,-1) 



respectively. (Normalisation is omitted for clarity.) 
Hence, the ground state Eq is totally bonding, (an A\ 
irreducible representation of C\ v ,) but has more weight 
upon the dimerised In2-Iri3 bond, thus corresponding to 
a dimer bonded to the single site Ini . This is a very dif- 
ferent result to that obtained 6,7 '— for the free Vp, where 
the lowest lying state was totally symmetric with equal 
weight on all sites. This new result for [Znj n -Vp] in- 
dicates that a dimer (or antidimer) may form even with 
only one electron present and no degeneracy. (The Jahn- 
Teller theorum does not state that symmetry cannot be 
reduced in the absence of degeneracy. ) This possible sym- 
metry lowering is weak, however: the energy gained is 
only ~ S 2 , the same order as the elastic energy lost. It 
comes about because the dimerisation/antidimerisation 
breaks the E irreducible representation of C^ v up into a 
B\ and a further A\ represention of C\ v . This new A\ 
representation can mix with the other ones, leading to a 
weak gain in energy, even though Jahn- Teller does not ac- 
tually require it here. Taking elastic energy into account, 
the total energy of the dimer structure for Q = +1 (with 
one electron) is 

E Q=+1 = -t(l + 5 + F{6)-3\5 2 ) (11) 

Taking the derivative of E® =+1 with respect to <5 we find 
that the lowest energy solution has S > (is dimerised) 
as long as < A < 0.17. This appears to correspond to 
a weak value of the elastic constant. For Q = — 1 — > —3 
a stable dimer state exists for any value of A and has a 
lower energy than the symmetric structure. Our DFT 
calcualtions indicate that the case of [Znj n -Vp] in InP 
is near to this critical value of A = 0.17: if we do the 
relaxation using a 4x4x4 k-point grid we do not find a 
stable dimer or a stable antidimer. However, if we use the 
less accurate 2x2x2 k-point grid we do indeed find a non- 
symmetric structural groundstate for Q = +1, namely a 
2.1% dimerised structure, with the symmetric structure 
an excitation 0.002 eV above. For Q = +0 we never find 
stable dimerised or antidimerised structures. 

Turning to the model excitations for the dimer struc- 
ture, E\ also corresponds to an In2-In3 dimer, (another 
Ai representation,) but now anti-bonded to the single 
site, upon which most of the weight is placed. E% is an 
antibonded dimer (an antidimer - a B\ irreducible rep- 
resentation of Ci v ) with no weight at all on the single 
site. For the antidimer structure (6 < 0) the groundstate 
remains the same (at least for small i5) but there is less 
weight on the dimer bond and more on the single site. 
Also, the order of the excitations is reversed, with the an- 
tidimer state (E3 above) becoming the lower of the two. 
If we now assume a form for 4>i related to sp 3 hybrids we 
can see what sort of charge densities we expect to see for 
the $j , as sketched in figure |H| 

In the case of the excitations for the symmetric struc- 
ture, we find 

*i, 2 = *± = (l,e ± ^,e ±4 f l ) (12) 
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with the two being, of course, degenerate. Unfortunately, 
describing the lifting of this degeneracy as a result of 
Coulomb and elastic effects (the two will almost certainly 
combine) lies well outside the scope of this model or pa- 
per. Hence what is shown in the figure is just the two al- 
ternatives with C3 symmetry and three node-like planes, 
based upon the form of $±. The one with low charge 
density along the In-In lines is likely to have a smaller 
overlap with - and hence Coulomb repulsion from - the 
symmetric E groundstatc. Hence it has been drawn as 
the lower of the two excitations. This cannot be pre- 
dicted with more certainty upon the basis of the current 
model calculation, however. 



TABLE VI: Model predictions for the lowest lying structures 
(marked L.L.) of [Znj n -Vp]^, together with structural exci- 
tation energies to leading order in S, for occupancies of to 
6 electrons (# e~). DM = dimerised, SY = symmetric, AD 
= antidimerised. 



# e~ Q DM SY AD 

+2 - L.L. 

1 +1 L.L. «f ^ 

2 L.L. 

3 -1 ^ St L.L. 

4 -2 Y 25t LX - 

5 -3 %i St L.L. 

6 -4 - L.L. 



With the obvious exception of the excitations for the 
symmetric structure, the similarity between the model 
charge densities in figure |5] and the ab initio charge den- 
sities in figure El is obvious. We are indeed seeing the 
predicted combinations of one electron orbitals from the 
three dangling In dangling bonds, controlled by Jahn- 
Teller effects which the model can also predict. 



Using equationElwe can now evaluate the energy of the 
symmetric, dimerised and antidimerised states for each 
value of electron filling and compare them. For exam- 
ple, with five electrons present ([Znj n -Vp] ) the sym- 
metric state has energy — t, the dimerised state has en- 
ergy — t(l + (5) and the antidimerised state has energy 
-f (1 + S + F(S)) « (-t(l +8+£). Hence the lowest 
is the antidimer, with the dimer as an excitation only 
~ higher, and the symmetric state lies ~ St higher 
still. Comparing to table El we indeed see that the an- 
tidimer is the most stable, followed by the dimer, with the 
symmetric structure the least stable. Table IvTl shows the 
leading terms in these excitation energies for fillings of 
to 6 electrons. Since both elastic and Coulomb energies 
favour the symmetric structure we anticipate only the 
symmetric structure for Q = +2 and —4, as observed. 
We also anticipate that the excitation energies for the 
symmetric forms will be lower than those indicated in 
the table. We find that for 1 or 2 electrons the dimerised 
structure to be the most stable, closely followed by the 
antidimerised structure. For 3, 4 or 5 electrons we ex- 
pect the antidimerised structure to be the lowest, with 
the symmetric structure always the least. The agree- 
ment with the ab initio results is very close - the only 
real exception is that for Q = — 3 our ab initio calcula- 
tions found the dimerised state to be slightly more stable 
than the antidimerised. The ab initio energy difference 
was tiny, however, over an order of magnitude less than 
the accuracy anticipated from such a calculation. On the 
other hand, the model does predict that the symmetric 
structure is at its least stable for charge Q = — 2 and this 
is indeed the only charge state for which we are unable 
to find the symmetric structure using DFT. 



VI. CONCLUSIONS. 

We have calculated the structure and transfer levels 
of the [Znj n -Vp] complex using ab initio DFT methods 
and have presented detailed results for the dependence of 
the structure and structural excitations upon the charge 
state of the complex and hence upon the Fermi level. We 
found that, for all charge states below -5, the Zn atom 
lies in a DX position, sp 2 hybridised, with an empty 2p z 
orbital directed towards the In triangle and very con- 
stant bond lengths to the surrounding P atoms. Elec- 
trons added (below Q = —5) go into localised orbtials 
built from the three In dangling bonds - thus reducing 
the In-In distances with each added electron. We also 
found that there is a simultaneous reduction in the Zn-In 
distances, which is of electrostatic in origin. This was due 
to the increased charge on the In triangle and the weaker 
screening of the Zn core along (111) owing to the empty 
Zn 2p z orbital. The combination of these effects leads to 
a systematic reduction in the effective vacancy volume 
of the complex for charge states from Q = +2 — * —4. 
This is in spite of the fact that the DX structure of the 
complex remains stable with rising Fermi level right up 
until a -4/-6 negative U transfer level, about 0.5 eV into 
the conduction band. For charge states Q = — 1 — * —3 
we also found complex Jahn- Teller behaviour, with up to 
three different stable structures - symmetric, dimerised 
and antidimerised - which are all very close to degen- 
eracy. For Q = — 1 the most stable structure (just) is 
dimerised, whilst it is antidimerised for Q = —2 and -3. 

In addition, we have presented a simple tight bind- 
ing model which is able to predict and explain all of the 
structural properties of the complex - including the ex- 
istance of three possible nearly degenerate structures for 
Q = — 1 — > — 3. Such a near degeneracy between dimer 
and antidimer structures has not been seen for the free 
Vp and this too the model can predict: the same model 
applied to Vp 1 , for example, shows a separation between 
dimerised and antidimerised structures on the order of 
5, compared to S 2 here for [Znj n -Vp] . (Details of these 
calculations will be presented elsewhere^) Success with 
such a simple model is very pleasing and may in the fu- 
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ture even help to narrow the range of initial configura- 
tions required when searching for structural minima in 
similar ab initio calculations. 

Comparing our results to those obtained experimen- 
tally using positron annihilation, we, as they, find that 
the most stable charge state of the [Znj n -Vp] complex 
is Q = in p-type material. We find a binding energy of 
0.39 eV which is also in the range they anticipate. Re- 
garding the structure of the complex, we, as they, find 
a larger effective vacancy volume than that of the free 
Vp. Upon excitation to the -1 charge state we, as they, 
see a reduction in volume. At 0.50 eV our 0/-1 transfer 
level is a little high for agreement, but, as was noted, 
these transfer levels are not reliable since the supercell is 
too small. Further studies of the complex in larger cells 
are required to check how this effects the position of the 
transfer levels, but our results seem otherwise to be in 
very good agreement with the experiments. 

Regarding the issue of compensation, our results con- 
firm that, should free Vp somehow be present in heav- 
ily Zn-doped InP, they would indeed form stable com- 
plexes with Zn. However, the formation energy we find 
for Vp is 1.96 eV, which is still rather large. Even for 
a phosphorus atom neighbouring a Zn dopant the effec- 
tive formation enery is 1.96 - 0.39 = 1.57 eV, which is 
not small. The positron annihilation experiments clearly 
demonstrate that at least some Vp are present in heavily 
Zn doped Czochralski grown InP. However, the lack of a 
low formation energy leaves it unclear whether, and per- 
haps unlikely that, Vp generation is the primary cause of 
compensation in InP:Zn. On the other hand, our DFT 
results indicate that the formation of [Znj n -Vp] com- 
plexes from existing Znj n and Vp is certainly not re- 
sponsible for compensation: starting with fully ionised 
Zn, the complex forming reaction is 

Zn In+ V P^[ Zn In- V p] ( 13 ) 

in which no free electrons or holes are created or ab- 
sorbed. The compensation issue here has to do with the 
formation of the phosphorus vacancies in the first place, 
and is not affected by the formation of the [Znj n -Vp] 
complex. Fully understanding the issue of carrier com- 
pensation in InP:Zn from an ab initio point of view may 
thus require more work. 



VII. FIGURE CAPTIONS. 

A. Figure 1. 

Structure of the supercell containing the neutral 
[Znj n -Vp] complex. Atom types are dark grey: In, pale 
grey: P, black: Zn. 



B. Figure 2. 

ELF plots for the neutral [Znj n -Vp] complex. For 
clarity most atoms are omitted and the edges of the Zn- 
In tetrahedron are drawn. Note, however, that there are 
no Zn-In bonds, a) ELF in a plane containing Zn, 1 In 
and the (111) axis, b) ELF in a plane containing Zn and 
2 In. 

C. Figure 3. 

Charge density difference plots for the [Znj n -Vp] com- 
plex in the Q = relaxed geometry, but with an ex- 
tra electron added (upper row) or removed (lower row). 
Shown are isocharge surfaces for the charge density dif- 
ference itself. The surfaces shown are for 80%, 50% and 
20% of the peak value. 

D. Figure 4. 

50% isosurface plots for charge density differences 
for adding (top row) and removing (bottom row) elec- 
trons from thcQ = +1, and -1 relaxed structures for 
[Zn In -Vp]. 

E. Figure 5. 

Top views (down the (111) axis) of the 50% isosurface 
plots for charge density differences corresponding to the 
transfer levels of [Znj n -Vp] from Q = +2 to -4, for the 
symmetric, dimerised and antidimcriscd structures. 

F. Figure 6. 

Side views of charge density difference plots for a) the 
Zn 2p z related orbital filled at the -4/-5 and -5/-6 trans- 
fer levels, b) the delocalised orbital emptied in removing 
an electron from Q = +5 and c) the orbital emptied 
in removing an electron from the mctastable shrunken 
structure for Q = +2. 

G. Figure 7. 

Irreducible representation splittings as symmetry is 
lifted from (9 3 to C\ v or C3. For C\ v the dimer case 
is shown. For the antidimcr the B\ level lies below the 
second A\ level. 

H. Figure 8. 

Sketch of the wavefunctions anticipated for the three 
structures of [Znj n -Vp] , based upon the model. (DM = 
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dimerised, SY = symmetric, AD = antidimerised.) 
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